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The electromagnetic field of the ultrarelativistic Reissner-Nordstr0m 
solution shows the physically highly unsatisfactory property of a van- 
. ishing field tensor but a nonzero, i.e., S- like, energy density. The aim 

' of this work is to analyze this situation from a mathematical point of 

view, using the framework of Colombeau's theory of nonlinear general- 
■ ized functions. It is shown that the physically unsatisfactory situation 

is mathematically perfectly defined and that one cannot avoid such sit- 
os : uations when dealing with distributional valued field tensors. PACS 

numbers: 04.40. Nr, 02.30. Sa, 04.90+e, 02.90. +p 

\Q '■ 

o : 

\o . 

On ■ 

^ , i. introduction 

I 

^Jy Recently, there has been some interest in gravitational shock wave geometries 

produced by ultrarelativistic particles Bui In field and string theory,qae is interested 
especially in high energy scattering processes in these geometriesJaQ 
k> | In 1971 Aichelburg and Sexfl derived the ultrarelativistic limit of the Schwarz- 

schild metric, which turns out to be a pp-wave with a distributional, i.e., <5-like 
profile fufi&tion. This result has been reproduced by several authors using various 
methodspEl but all of them invoking distributional techniques. Lousto and SanchezO 
derived the ultrarelativistic limit of the Reissner-Nordstr0m solution. To gain a 
distributionally well defined limit of the metric they had to rescale the charge by 
e 2 = — w 2 ) 1 / 2 which also forces the limit of the electromagentic field to vanish 
(even in the space of distributions). However the limit of the energy stress tensor of 
the eletromagnetic field is nonvanishing, thus leaving us with the physically highly 
unsatisfactory situation of a vanishing field producing a (5-like energy density. 

The aim of this work is to discuss this situation from a more mathematical point 
of view. Computing the electromagnetic stress tensor from the field tensor is a non- 
linear operation and thus not defined within the framework of classical distribution 
theory. Thus, strictly speaking, if one is dealing with distributional-valued ficM 
tensors one has to go beyond distribution theory. In the 1980s J.F. ColombeautiiI9 
developed a theory of generalized functions providing the possibility of a prod- 
uct of distributions. He constructed differential algebras Q containing the space 
of smooth functions as a subalgebra and the space of distributions as a subspace. 
Thus Colombeau's theory provides us with a natural framework for discrj^sing situ- 
ations like the above-mentioned one from a more abstract viewpoint Jld'E3 Moreover 
most recently first applications o£,this very framework to problems of classical field 
theoriesEJ and general relativityE.il have appeared. 
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This work is organized in the following way. In Sec. II we briefly recall the basics 
of Colombeau's nonlinear theory of generalized functions. To get some working 
knowledge, in Sec. Ill we take a look at the ultrarelativistic limit of the Coulomb 
field in flat space, and finally in Sec. IV we discuss in some detail the electromagnetic 
field of the ultrarelativistic Reissner-Nordstr0m solution and its energy stress tensor. 

II. BASICS OF COLOMBEAU THEORY 

Let f2 be an open subset of R™ or a smooth manifold. We denote by V(il) the 
space of test functions on fi, i.e., the space of smooth functions with compact 
support in f2, and by X>'(f2) the space of distributions on Q. Finally the action of a 
distribution w on a test fuction ip we denote by (w,ip). 

In V no meaningful product can be defined. Moreover, as L. Schwartz showed in 
1954, there even exists no associative and commutative differential algebra contain- 
ing the space of continuous functions as a subalgebtartiiat allows a linear embedding 
of the space of distributions in it. J.F. ColombeamV^i introduced differential alge- 
bras Q containing the space of distributions as a subspace, and the space of smooth 
functions as a faithful subalgebra, thus providing a natural framework for studying 
nonlinear operations on singular data, i.e., distributions. 

In this paper we consider the "special variant" of Colombeau's algebra, whose 
definition we briefly recall. Set 

£ M (Sl) = {«)e G (C°°(n)Y ^ : VK CC nVa G N%3N > : 

sup | d a u e (x) \— 0(e~ N ) (e->0)}, (1) 

jV(fi) = { {u e ) e E (C°°(n))^ : VK CC n\/a E N£,VM > : 

sup I d a u e (x) |= 0(e M ) (e-fO)}. (2) 

xeK 

£m(Q) is a differential algebra with pointwise operations and Af(fl) is an ideal in 
it. We define the algebra of generalized functions, or Colombeau algebra, by 

g(n) := e M (n)/Af(n) (3) 

and denote its elements by 

U = (u e )e + A/"(fi) . 

Distributions with compact support on 1" can now be embedded into C?(R ra ) by 
convolution with a mollifier p e , defined as follows; let p € 5(R") (Schwartz's space) 
with the properties J p(x) dx = 1 and J x a p(x) dx = Va 6 N™, | a \> 1, then we 
set p e (x) := (l/e n )p(x/e). So we have the map l (uj) = (uj * p e ) e +ftf(M. n ). 
This embedding can be "lifted" to an embedding i : X>'(f2) ^ G(&) by means 
of sheaf theory while smooth functions are embedded as constant sequences, i. e. 

*(/) = (/)*• ' 

Next we briefly recall the concept of association in the Colombeau algebra, of 
which we make extensive use in the physical calculations of Sees. Ill and IV. One par- 
ticular property of Colombeau's algebra is the so-called "coupled-calculus," namely 
equality in the algebra of generalized functions and equality on distribution level. 
The latter is defined as follows: a generalized function u is called associated to 
(u w 0) if for one (and hence any) representative (u e ) e we have 

lim / u e (x)ip(x) dx = 0, Vip € T>(Sl). (4) 
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Association defines an equivalence relation in the Colombeau algebra that clearly 
is coarser than equality and is compatible with differentiation, i. e., u rs v => d a u « 
d a v but is not compatible with the algebra multiplication, i. e., u « v ^> wu sa wv. 

If a generalized function u is associated to l(w), where w is a distribution, then 
one says that w is the distributional shadow or the associated distribution of u. Not 
all Colombeau functions have a shadow-for example S 2 , i.e., t(<5) 2 -but if a generalized 
function has a shadow the latter is unique. Different generalized functions may have 
the same shadow, providing us with the possibility to model singular situations in 
various waysJlj 



III. THE ULTRARELATIVISTIC LIMIT OF THE COULOMB FIELD 

In Cartesian coordinates the Coulomb field is given by the Maxwell field tensor 
components F 0a = — (ex a )/r 3 where r 2 = x 2 + y 2 + z 2 and a — 1,2,3, and the 
other components vanishing. We apply a Lorentz boost with velocity v in the x 
direction and, in order to simplify computations, we additionally transform to a 
frame associated with null coordinates; in particular 

t i — ► i = j(t + vx) i— > u = x — i , 

x i-^ x — j(x + vt) i— > w = x + i , (5) 
where 7 := (1 — v 2 )~ x / 2 . Thus the boosted field tensor takes the form 



F 



ik 



e(l 



[[x - vt) 2 + (1 - v 2 )p 2 } 3 / 2 



I 2(x-vt) (l-v)y (l-v)z 

-2(x-vt) -(l + v)y -(l+v)z 

-(l-v)y (l + v)y 

V -(l-v)z (l + v)z 



where p 2 = y 2 + z 2 . 

We could now view the components of the field tensor as sequences of locally inte- 
grable functions parametrized by the boost velocity v, and calculate their limits in 
the space of distributions as v — > 1. However, it is well known that in this limit some 
components of the field tensor tend to 8{u) and hence the energy momentum tensor 
Tik = FuF l k — \/AF m iF lm of the ultrarelativistic field cannot be calculated directly 
due to the lack of a multiplication in the space of distributions. Alternatively one 
could think of first computing the energy momentum tensor of the boosted field and 
then taking the limit. However, in this approach some components of the momen- 
tum tensor diverge in the space of distributions, thus making it again impossible to 
compute the energy momentum tensor of the ultrarelativistic Coulomb field within 
the framework of classical distribution theory. 

To overcome these difficulties we embed the components (F lk v ) v of the field 
tensor into the Colombeau algebra. The physically most natural way to do so 
is to view a whole sequence (for example (F 01 „)„) as one Colombeau function. 
Generalized functions in the sense of Colombeau are classes of sequences of smooth 
functions, but our sequences are not even defined on the set A = M. 4 \ {x = vt (0 < 
v < 1) A p = 0}, where from now on we neglect the bars over x and t to simplify 
the notation. So in order to get an "embedding" of (F lk v ) v into Colombeau's 
algebra without cutting the whole set A from the domain we employ the following 
construction that only forces us to restrict the domain to ft := R 4 \ {x = t A p = 
0}, which is physically reasonable since we precisely cut out the world line of the 
ultrarelativistic particle. Let X v £ C°°(R) (0 < v < 1) be such that 

x v (o = i, m>4(i-«) 2 , 

= 0, \t;\<(i-v) 2 . 
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Using the abbreviation a := (x + vt) 2 + p 2 we define 

F ik v = 




a > 4(1 - v) 2 
4(1 -v) 2 > a > (1 - v) 2 
(l-v) 2 >a, 



(6) 



which lies in Q(Vt) and equals F lk v outside a cylinder of radius 2(1 — v) with axis 
x = vt in the p = plane, whereas it is smoothed down to zero inside a cylinder 
with half this radius. Note however that F %k v is smooth even on M 4 . Now F lk v 
is the best possible approximation for F %k v in Q(fl) in the following precise sense: 
Given any compact K C f2, finally (i.e., v large enough) F lk v equals F lk v on an 
open neighborhood of K. 

Putting u(f) = 1 - e we get the desired "embedding" of the field tensor into the 
Colombeau algebra 

(F ik v )v - (F ik e) ) e +Af(n). 

Now we are able both to compute the energy momentum tensor in Q (SI) , and to 
make use of the concept of association, which corresponds physically to taking of 
the ultrarelativistic limit. First we calculate the association relations for the field 
tensor. 



Proposition 1 



/0 



y 

\0 z 





—y —z 



4e 

5^) 



(7) 



Proof: We have to compute the following limits 



lim / F l v k e) (t,x,y, z) ip{t,x,y, z) dtdxdydz , 



(8) 



where ip is an arbitrary test function in 2?(il). We only carry out the calculations 
for F 01 and F 21 as the other components are either related to these by symmetry 
or the limits can be taken in a closely analogous way. We choose v so large that on 
an open neighborhood of the support of tp we have F* k — F lk v , and for simplicity 
we set [ ] := (x — vt) 2 + (1 — v 2 )p 2 . Then we have 



^.01 



suppM 



<2e|M| c 



(1 - v 2 ) (x - vt) 



13/2 



dx 



dtdydz 



— 8ne\/l — i^H^Hoc J J - pdpdt 



(9) 



which vanishes in the limit, since the remaining integral has to be taken over a 
compact set only. 

To prove F 21 v ^ « (Aey/p 2 )5(u) we first note that for f v ( e ), defined from /„ = 
ey(l + v)(x - vt)/(p 2 [ ]V2) in the same way as F %k v ^ is defined from F lk v we 
have 



lim 



supp( v ) 



2 jy[e(x-t)-6{t-x)],y 



(10) 
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where we have used Lebesgue's dominated convergence. Differentiating this relation 
with respect to x we get the desired result. □ 

Next we compute the energy stress tensor in Q{Vt) by componentwise multiplica- 
tion to get 

T ik v(e) 



2((x - vt) 2 + (1 - t) 2 )p 2 ) 3 

/ ±(l + t>)V \{x-vt) 2 (l + v)(x-vt)y (l + v)(x-vt)z \ 

±(x-vt) 2 ±(l-v) 2 p 2 -(1 - v)(x - vt)y -(1 - v)(x - vt)z 

(l + v)(x-vt)y -(1 - v)(x - vt)y (x - vt) 2 + -2(1 - v 2 )yz 



. (1 + v)(x - vi)z -(1 - v)(x - vi)z -2(1 - v 2 )yz (x-vif + 

\ +(l - v 2 ){y 2 - z 2 ) ) 



(11) 



for a > 4(1 — v) 2 and smoothed down to zero "inside" in the same way as in the 
case of F lk v ^y 

Proposition 2 

f 4fc «0 V(z,fc)^(0,0) 

T o has no associated distribution. (12) 



Proof: The calculation for all the components is closely analogous to the one given 
above with the exception of Too, so we shall only examine this one in detail. As a 
test function ip we take 1/p 2 on the set x,fe [— N, N], 4> e [0,2n] and p e [?"o,ri] 
(0 < r ), where <p is the polar angle in the (y, 2;)-plane, and let it smoothly approach 
zero "outside" . Thus we have for large enough v 



I 



N N 2w i-i 

e 2 (l + v) 2 (1 - v 2 ) 2 f f f f p 2 1 
T QO v(e)ip> j / / / / r~13~ pdpd<j)dxdt 



SUpp( V ) -AT -N r 

\n,N,N 



e 2 7r(l + v) 2 (x — vt) x — vt 



16 vp*Vl^ arCtan VT^P ]r °>- N >- N 

e 2 7T(l+v) 2 ])]ri ,N,N 

32 V P 2 U l)l r ,-N,-N 

e 2 7r(l +v) 2 ( (x - vt) 2 \ NiN 



32 vp 2 ^ \ X + (1 - V 2 )p 2 ) \ro\-N,-N ■ (13) 

In the limit the second term diverges whereas the other two terms vanish. □ 

We see that computing the energy momentum tensor we fall out of the class of 
Colombeau functions with associated distributions. Taking the square of the two 
functions {F^ 2 ^) e and (F^ e ^) ei both of them associated to the S distribution, we 

get the generalized function (T 00 „ (e) ) e = -(1/4) {[F^) 2 + (T^) 2 ) not allowing 

any associated distribution, showing that (T o„( e )) e is an object only defined in 
the Colombeau algebra. Due to the properties of Colombeau's theory one could 
subject it to a wide class of nonlinear operations and even use it as a source term 
for nonlinear PDEs and then check if the resulting generalized function has an 
associated distribution. 
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Of course this result mirrors the fact that if we had computed the energy tensor 
"naively" from the ultrarelativistic limit of the field tensor in V we would have 
obtained T o — (4e 2 / p) u S 2 (u)," all other components vanishing. This suggests a 
"commutativity" of the limiting procedure and "multiplication," i.e., the compu- 
tation of the stress tensor from the field tensor. However, we again point out that 
no meaningful product can be defined within the framework of distribution theory 
but only in the Colombeau algebra. As a limiting procedure we have used the con- 
cept of association in G(Jl), which is not compatible with the product, and so the 
above-mentioned "commutativity" does not hold in general, as we are going to see 
in the next section. 



IV. THE ULTRARELATIVISTIC LIMIT OF THE 
REISSNER-NORDSTR0M SOLUTION 

The RN metric representing the gravitational field of a point particle with mass 
m and charge e in isotropic coordinates is given by 

ds 2 = (l-^- + -4^) dt 2 - (1 + ™ + ^^) 2 + d ^) > ( 14 ) 
v r(r) r 2 (r)' v r 4r 2 ; v ; v ; 

where r = f[l + m/f + (m 2 — e 2 )/(4f 2 )] is the radial coordinate in a Schwarzschild- 
like coordinate system. We apply a boost in x direction relative to an asymptotic, 
Cartesian coordinate system associated with the isotropic radius f. To obtain a 
distributionally well defined result in the ultrarelativistic limit we have to rescale 
the mass and the charged in the following manner: 

m = p/j, e 2 = p 2 /7, 

where in the limit, we keep the momenta p and p e fixed. The rescaling of the mass 
is well motivatedQ and saves the energy of the particle from diverging due to its 
finite rest mass by keeping the total energy p fixed and letting m approach zero 
in the ultrarelativistic limit. The rescaling of the charge is not too well motivated. 
However, the limit of the metric is given byu 

ds 2 = \ 8p\np H — 7T— > 5(u) du 2 — dudw ~ dy 2 ~ dz 2 , (15) 
I 2 p 1 

where u and w again denote null coordinates defined as in the last Section. This 
metric belongs to the class of pp waves and is flat everywhere except on the null 
plane u = which is normal to the boost direction and contains the particle. 

However, the rescaling of the charge dramatically acts on the limit of the elec- 
tromagentical field and its stress energy tensor. In a frame associated with null 
coordinates the boosted field tensor takes the form 

(0 x — vt y z \ 

-x + vt -vy -vz I 

-y vy o o ' 

-z vz / 

where the factor G(f) — 4f 2 (4f 2 — m 2 + e 2 )/(4f 2 + Arm + m 2 — e 2 ) 2 arises from the 
transformation to isotropic coordinates and we again have neglected the bars over 
x and t. Now we embed (Fik v ) v into the Colombeau algebra in the same manner as 
in the Coulomb case, again denoting the G(Q) function best approximating by 
(F ik u( e ))e- By calculations very similar to those carried out in Sec. II and taking into 
account the rescaling of mass and charge we get the following association relations. 
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Proposition 3 

F lk « Vi.fc. (16) 



Next we compute the stress energy tensor. For a>4(l — w) 2 we have 
- P 2 G 2 (f) 

4fc ((a; - vi) 2 + (1 - v 2 )p 2 ) 3 1 U 



/ (1 - u) 3 / 2 (l + t)) 7/2 ( 3 11 - s 00 )(x - i,t) 2 /4 (1 - v 2 f' 2 (g aa + S ")(x - trf) 2 /4 

+(1 - „)«/a ( i + „)9/ 2 ( s °° + fl ")p"/4 +(l - « 2 ) 7 /2 (s ii _ 9 oo ))9 2 /4 

(1 - „ 2 )5/ 2 ( s °° + s ")(x - «t) 2 /4 (1 - «) 7 / 2 (l + „) 3 / 2 ( S 11 - g°°)(z - ^) 2 /4 
+ (1 - « 2 ) 7 / 2 (g 11 - S 00 )p 2 /4 +(1 - „)»/*(! + vf'H.g oa + a 11 )^^ 

(1 - vf' 2 {l + v) 7 ' 2 g aa (x - vt)y -(1 - v) 7 ' 2 (l + vf' 2 g aa (x - vt)y 

(1 - vf' 2 {l + v) 7 ' 2 g aa (x - vi)z -(1 - v) 7 ' 2 {l + v) 5/2 g 00 (x - vt)z 



{l-vf' 2 (l + v) 7 ' 2 g aa {x-vt)y (1 - vf' 2 {l + v) 7 / 2 g°°(x - vt)z \ 

-(1 - t-) 7 / 2 (l + vf' 2 g 00 (x - vt)y -(1 - «) 7 / 2 (l + vf' 2 g 00 (x - vt)z 

(1 - v 2 f' 2 g Q0 [(x - vt) 2 + (1 - v 2 )(z 2 - y 2 )] -2(1 - v 2 ) 7 ' 2 g 00 yz 

-2(1 - i, 2 ) 7 / W (1 - « 2 ) 5/ V°[(* - vt) 2 + (1 - t, 2 )(y 2 - Z 2 )] / 



where 

-.2 _ „2 c 2 _ ^,2 

\2/ 



5 0U w = a + T + ^-2-ra + -«5-)- J ( 18 ) 



and 



S 11 (r) = (l + f + ^^)- 2 (19) 

and Tifc smoothed down to zero "inside" again in the same way. By a rather lengthy 
calculation, expanding G(f), g oa and g 11 and again taking into account the rescaling 
we get the following. 



Proposition 4 



T lfc «0, V(i,fc)^(0,0), 



V. CONCLUSION 

The fact that the "square" Too of the generalized functions (i r ifc) e associated to 
is not associated to but to the S distribution is not surprising from the viewpoint 
of Colombeau's theory. Association is not compatible with multiplication in the 
algebra of generalized functions. Thus we can say that the physically unsatisfactory 
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situation of a vanishing field but nonzero stress tensor is mathematically perfectly 
defined. We sec the incompatibility of linear idealisations like the S distribution and 
nonlinear computations. Colombeau's theory clearly cannot solve this principal 
defect but provides us with a language and framework that makes it possible to 
discuss such situations. 

The question of the physical relevance of a zero field producing a nonvanishing 
energy density clearly has to remain open. But from the mathematical viewpoint 
such a situation cannot be avoided if one is dealing with distributionally shaped 
fields. 
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